arXiv:1508.04656vl [astro-ph.SR] 19 Aug 2015 


Title : will be set by the publisher 
Editors : will be set by the publisher 
EAS Publications Series, Vol. ?, 2015 


MASS LOSS OF RED SUPERGIANTS: A KEY INGREDIENT 
FOR THE FINAL EVOLUTION OF MASSIVE STARS 

Cyril Georgy 1 and Sylvia Ekstrom 2 


Abstract. Mass-loss rates during the red supergiant phase are very 
poorly constrained from an observational or theoretical point of view. 
However, they can be very high, and make a massive star lose a lot of 
mass during this phase, influencing considerably the final evolution of 
the star: will it end as a red supergiant? Will it evolve bluewards by 
removing its hydrogen-rich envelope? In this paper, we briefly sum¬ 
marise the effects of this mass loss and of the related uncertainties, 
particularly on the population of blue supergiant stars. 


1 Introduction 


Mass loss is an important driver of the evolution of massive stars. Stars more 
massive than ~ 15 Mq lose du ring their life, ~ 4 — 7 (for a 15 M^) to ~ 
100 Mq (for a 120 Mq model, see lEkstrom et ah . 201 2t Ictiieffi and Limongi . l2013l) . 
This loss of matter will progressively uncover deep layers of the star, modifying 
the surface chemical composition together with internal mixing processes. Stellar 
evolution models are currently not able to compute the mass-loss r ate consistent ly, 


and thus rely on recipes, both observational or theoretical (e.g. Ide Jager et al 
Il988l: Kudritzki and Pula . 2000t Vink et al. . 2000h 


One of the worse constrained of these rates is the red supergiant (RSG) one. As 
shown in Fig. |I] the mass-loss rates have a huge scatter (three orders of magnitude). 

In the latest release of stellar models grids published b y the Geneva stellar 
evolution group ( Ekstrom et all 2012t Georgy et all l2013lb stars more massive 
than about 20 Mq have layers in their envelope that have luminosity considerably 
higher than the Eddington one during the RSG phase. This could ease the superfi¬ 
cial layers to escape the star. We have thus decided to increase the standard RSG 
mass-loss rates by a factor of 3 in our computations for stars more massive than 
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Fig. 1. Mass-loss rate during the RSG phase. Observations from iMauron and Josselinl 
ll201ll l (blue stars) and Ivan Loon et al.l (|2005l l: M-type stars (red stars), MS- and S- 
type stars (green), and carbon RSGs (purple). The black curves show the mass-loss 
rate according to Ide Jager et al.l l|l98^ 1. for log(T e fr) = 3.5 (solid), 3.6 (dashed) and 3.7 
(dotted-dashed). The red zone correspond to the mass-loss rate used in the Geneva stellar 
evolution code for standard mass-loss rate (lower line), or increased one (top line). 


20 M e . Th is is shown by the red diagonal band in Fig. [I] where the classical mass- 
loss rate of lde Jager et al. ( 198 8) is also plotted for comparison for three different 
effective temperatures. Our mass-loss rates are well in the range of the observed 
ones. In the following, we will discuss a few conseq uen ces o f such an increased 


mass-loss rates during the RSG phase (see also Vanbeveren et al. . 1998h . 


2 Second crossing of the HRD 

As a direct consequence of the high mass-loss rates during the RSG phase, the 
stellar models lose much more easily their extern al h ydrog en-rich envelope. This 
favours a bluewards motion in the HRD ( Giannonei . Il967 1. In our models, this 
occurs for star more massive t han 25 Mq (non-rotating models) or 20 Mq (rotating 


luminosity of Galactic RSGs ( 

Ekstrom et al. 

2012), as well as the maximal mass 

of observed type IIP supernova progenitors 

Smartt et al., 2009j). 


A first consequence of this second crossing of the HRD is the co-existence of 
two populations of blue supergiants (BSGs), one from stars just after the main 
sequence, going to cross for the first time (group 1), and the other one from stars 
coming back from the red side of the HRD (group 2). The study of the pulsational 
properties of both groups shows t hat they are very different, and can thus be 
used to determine the membership (iSaio et all 1201311 . Figure [2] shows the excited 
periods that are visible at the surface. Before the RSG phase (left panel), group 
1 stars have only a few excited modes and almost none in the log(T e ff) range from 
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Fig. 2. Left and central panel: HRD (on top) and excited pulsation periods visible at 
the surface of our models of the group 1 (left) and group 2 (centre). Observed positions 
of some a Cyg variables are indicated. Adapted from ISaio et al.l l2013l l. Right panel: 
Position of observed BSGs in the FWL relation plot, and expected positions of a synthetic 
population of rotating stars (coloured region). Adapted from iMevnet et al] (l2015bl') . 


4.0 to 4.2. On the other hand, group 2 stars (centre panel), that have lost a huge 
amount of mass during the RSG phase, have a much favourable L/M ratio, that 
allows for much more modes to be excited. The predicted periods are in good 
agreement with the observed period of a Cygni variables. This group of BSGs can 
also b e used to test the physics of convection in massive stars. We have shown in 
Georgy et al. ( 2014 ) that using the Ledoux criterion instead of the Schwarzschild 
one provides surface abundances for a Cyg variables in better agreement with 
ob servations. _ _ 

iKudritzki et al. ( 20081 ) have shown that observed BSGs lie in a linear relation 
if they are plotted in the plane log(L) vs. log(g/T e 4 ff ) where g is the surface gravity 
(see the stars in the right panel of Fig. [2]). Th is relatio n is called the flux weighted 
gravity-luminosity relation (FWLR). In IMevnet et al. ( 2015bl) . we have shown 
that our models are able to reproduce reasonably well this relation (right panel of 
Fig. [2]) - However, we may generate too much group 2 stars, producing a slightly 
off-relation region where we expect from our models a high density of BSGs that 
are not observed. This could indicate that our models have an initial velocity that 
is too high compared to the average velocity of a typical massive star, or that our 
enhancement of the mass-loss rate during the RSG phase is too high. 


3 Conclusion 


The two examples discussed above show the consequences of the mass-loss rates 
during the RSG phase on the further evolution of the star. Observations of BSGs 
provide constraints on the mass-loss rate that have to be used in the stellar evo¬ 


lution models during that phase (see also iGeorgvl 120121 ; Mevnet et al-, I2015al) . 


However, a better determination of the RSG mass-loss rates are extremely needed 
in order to improve the stellar modelling of the advanced stages of stellar life. 
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CfflAVASSA: Do you take into account the mass loss temporal variability during 
the RSG phase? As well as the chemical composition of the star? 

Georgy: The mass-loss rate recipe for the RSG phase depends usually on stellar 
parameters, basically the temperature and luminosity. So the long-term variability 
due to the secular evolution of the star is accounted for through that way. The 
possible short term variability is not accounted for. Depending on the typical vari¬ 
ation timescale, it could or could not be included in the stellar evolution codes: we 
need to have timesteps not too small. In the latter case, we would require a time- 
average mass-loss rates. Concerning the chemical composition, as far as I know, 
there is no mass-loss prescriptions that account explicitly for it in the literature. 

MORAVVEJI: We should not ignore the fact that to each datapoint on the FWLR 
contour plot, there should be attached to it a large error bar. So, in spite of that, 
is there a “statistically” significant difference between your predictions for rotating 
vs. non-rotating tracks? 

Georgy: First of all, I would like to precise that we have included an error in our 
contour plot, that is 0.15 dex in magnitude, and 0.075 dex in log(T e ff). We did not 
perform any statistical test to compare our predictions. They are slightly differ¬ 
ent, because the tracks are definitely different between rotating and non-rotating 
models. However, we are far to have a big enough observational sample to test the 
detail of the models so far. 


De Marco: Ledoux criterion is very different. It seems like we should have no¬ 
ticed that the Schwarzschild criterion was producing very different results. How 
likely is that this is the right explanation? 


Georgy: At the beginning, we used the classical Schwarzschild criterion. But 
quickly, we had to face the problem that even if our prediction about the pulsation 
peri ods were re l ativel y good, our predicted surface abundances were completely 
out ( Saio et all 120131) . So, we tried to change different things (overshoot, rota¬ 
tion, ...), but at the end, only the Ledoux criterion was able to reconcile both 
observations. However, and I want to be clear here, it does not mean that this 
is the correct solution. The treatment of convection in stellar codes is still very 
primitive, and it is well possible that the correct explanation is still different. 


Moffat: Last week some of us were at the Wolf-Rayet workshop in Potsdam. At 
that meeting it was noted that from Geneva models, WR stars come from stars 
with Mi n i > 25 Mq while RSG come from M- In [ < 20 M 0 , i.e. without overlap. 
Yet it was noted that the dense, rich cluster Westerlund 1 has a large population 
of both RSG and WR stars. What could be the reason to this? 


Georgy: It is true that according to our models, the overlap is very short in 
duration (but still exists). One important point to remember here is that a non- 
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negligible fraction of WR stars are probably produced in interacting systems. That 
could well be the case here. 



